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Introduction. The possibility of using cosmic ray induced neutrinos to 
detect oscillations in deep underground experiments has been considered 
by many authors. > Most of the estimates are made neglecting matter 
effects, or in the two neutrino case. The matter effects, however, are 
non-negligible : > in the two neutrino case, they reduce a mixing angle 
of 45° to 7.5° for 1 GeV neutrinos of squared mass difference 10~^ ev” 
going through the earth (average density = 5.52 g/cm 2 3 * ) making the oscil- 
lation totally unobservable. Also, they produce a natural oscillation 
length of about 6000 km in the case of massless neutrinos, invalidating 
all estimates based on wavelength arguments. Adding a third neutrino 
flavor considerably modifies the oscillation pattern and suggests that 
scales down to 5 x 10 -5 eV^ could be observable even when we take into 
account matter effects and the electron contribution to the incoming flux. 

In this paper, after defining the model of the earth used in the 
calculation, we will show the effect of matter on the probability curves 
for different cases, varying the masses and the mixing matrix. We will 
then give the prediction for the ratio upward V + v/downward V + V as a 
function of the zenith angle at Cleveland, neglecting angular smearing 
and energy threshold effects. 

1. Simple Model for the Earth: The 3 neutrino oscillation formalism 
has no analytic solution in the case of a continuous medium with varying 
density. 7 We therefore approximated the real density profile of the 
earth by a 3 layers model. The earth's density, especially in the core, 
is rather poorly known. 8 The only firm constraint, i.e. - that the 
average density is 5.52 g/cm^ and the diameter 12756 km, has been 
included in our model: the three layers have respective densities of 
3.5 g/cm 5 * , 5.5 g/cm 5 and 10.17 g/cm 5 and respective radial thicknesses of 
928 km, 1950 km and 3500 km. These densities and thicknesses come from a 
crude estimate made on one of the models quoted in Reference 6. 

2. Probabilities. We now study oscillations through a diameter of the 

earth using the notations of Reference 7. 

The main problem is the choice of a mixing matrix: in the case 

of Dirac neutrinos, the mixing matrix depends on 3 angles and one phase. 

We do not know of any model predicting a 3 neutrino mixing matrix with 
mixing angles large enough to be observed. 

Accordingly, we consider the following matrices, with U ai the ampli- 
tude to go from a charged-current eigenstate |v > to a mass eigenstate 
v. > of mass mi: a 

1 i 1 
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Case 1: 2 neutrino mixing matrix "Cahibbo like": 
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U 2 " 


cos a -sin a 
sin a cos a 




with a = 45° (maximal vacuum mixing) 

22.5° (lifting the Sm 2 (1) 

sign degeneracy' 7 ) 


Case 2: 3 neutrino mixing matrix with maximal muon mixing: 

This reduces the outgoing v u flux most. 
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Big 3-neutrino mixing with a phase. To study the phase-depend- 
use the Kobayashi-Maskawa matrix with a phase 6 and all angles 
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We study the probability P a to get flavor a from a beam of v e or 
when we let it evolve through matter and assume one of the previous mix- 
ings. The results are shown in Fig. 1 for 1 GeV neutrinos. 

We see that a squared mass difference of 10~4 eV 2 in the 2 neutrino 
case is unobservable due to matter effects. 3 > 4 However, this estimate 
is violated in the 3 neutrino case: for an appropriate mixing matrix, 

<5nr scales of order 10“ 4 eV 2 are attainable. Furthermore, graphs (b) 
and (c) clearly show that the oscillation pattern is totally modified by 
matter effects. Notice that the probabilities do not go back to their 
original value in the course of the oscillation. This is a matter effect 
which one always encounters in the presence of more than two uniform 
media. Also, we see that 3 neutrino oscillations are very sensitive to 
the phase of the mixing matrix and that matter effects can, in principle, 
distinguish between 6 wr > 0 and S m 2 < 0 as has already been pointed out. 


3. Upward to Downward Ratios at Cleveland. We now present results for 
the upward to downward ratio of contained events at Cleveland, corre- 
sponding to incoming fluxes partitioned into zenith angle bins such that 
.1 x n < cos 0 :£ .1 x (n+1) where 0 is the zenith angle and ne[0,9]. The 
energy is integrated for each bin from .2 GeV to 3 GeV. We neglect the 
effect of the layer above the detector and define the geometry assuming 
the detector at the surface of the earth. The zenith angle dependent 
fluxes that we use have been described elsewhere. 2,9,10 calculation 

is made at solar minimum and we assume a ratio v/v of 1.3. We take a 
fiducial volume of 2 x 10 33 nucleons* 1 and the inclusive cross sections 2 : 

= .737 x E x 10~ 3 ® cm 2 and o- =. E x .31 x 10“^ cm 2 with E the neutrino 
energy in GeV. We show in Fig.2 the results of our calculations in the 
form of ratios upward v + v /downward v + v. We first exhibit the 
effect of matter on 2 neutrino oscillations at the scale <Sm 2 = 10 -4 eV 2 . 
Although one would expect a big depletion (Fig. 2d) of the v flux from a 
naive argument based on vacuum oscillations, we see in Fig. 2c that mat- 
ter effects completely modify the picture and make the phenomenon unob- 
servable in this case. We also show results in the 3V case which sug- 
gest that it is possible to observe 6m 2 scales of order 10 -4 eV 2 . 
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4. Conclusions: We have shown in this work that matter effects totally 

modify the expectations for neutrino oscillations through the earth, 
making the accessible range for the squared mass difference of the order 
° f I°~ o eV t0 I 0 " 0 evZ in the 2 neutrino case and of the order of 
10" eV 2 to 10- 5 eV 2 in the 3 neutrino case. In both cases, the mixing 
needed to reach this precision must be unrealistically big. 

JwJ'k? energy averaging of § 3, we have taken the same electron and 
muon thresholds. An interesting possibility exists for the case where 
the electron threshold is lower than the muon threshold, as for the 
example of contained interactions shown in Fig. 2. One can then imagine 
a situation m which the upward v e interaction rate is increased without 
a significant decrease in the corresponding rate of Vy interactions. 

(This is possible because, at a given neutrino energy, v../v e ~ 2) . An 
oscillation effect maximum around 500 MeV could therefore increase the 
v e interaction rate via Vy •+■ v e more than it would be depleted via 
*£ v y wit ^out affecting the Vy interaction rate significantly since 
the latter is due primarily to neutrinos of energy E > 1 GeV. 
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FIGURE 2. 


FIGURE 1. Neutrino oscillations through a diameter of the earth. 
Horizontal axis: distance traveled through the earth; Vertical axis: 

probability P a . 

(a) Case 1 with | 6jn 2 | = lO - ^ eV 2 t The three curves give P e for a pure v 
incoming beam, for a = 45° (curve a), a = 22.5° and 6m 2 < 0 (curve b), 1J 
a = 22.5° and 5m 2 > 0 (curve c). 

(b) Case 3 with 6=0" m 2 = m 2 = 10" 8 eV 2 , m 2 = 10~ 3 eV 2 . 

(c) Same as (b) considering the earth as a vacuum. 

(d) Case 3 with 6 = 270°, m 2 = 0, m 2 = 1.5 x 10" 4 eV 2 , m 2 = 5 x 10" 5 eV 2 . 

(e) Same as (d) but with 6=0°. 

(£) Case 2 with m 2 = 0, m 2 = 10“ 4 eV 2 , m 2 = 5 x 10' 4 eV 2 . 

FIGURE 2. Ratios - , .CvlV^ u P . f or jo even bins in cos 6. 

(v+v)down 

Horizontal axis: cos 0; vertical axis: ratio . 

(a) Without oscillations. own 

(b) Case 2 with m 2 = 0, m 2 = 10 -4 eV 2 , m 2 = 5 x 10 -4 eV 2 . 

(c) Case 3 with <S = 270°, m 2 = 0, m^ = 1.5 x 10 -4 eV 2 , m 2 = 5 x 10~ 8 eV 2 . 

(d) Case 1 with a = 45°, 6m 2 = 10 -4 eV 2 , neglecting matter effects. 

(e) Same as (d) but with matter effects. 

(f) Same as (e) but with a = 22.5°. 










